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Chapter 1 

Adaptive Space-Time Beamforming in 
Radar Systems 



1.1 Introduction 



Space-time 

adaptive processing (STAP) techniques [Kleminl |2002h ]. |Melvin I |2004[ )] 



have been thoroughly investigated in the last decades as a key enabling 
technology for advanced airb orne radar applications fol lowing the seminal 



work by Brennan and Reed Brennan and Reed I (jl973h ]. A great deal of 



attention has been given to STAP algorithms and different strategies to 
design space- time beamfo r mers to mitigat e the e ffect of clutter and jam- 



ming signals [Reed et all (|1974t )]- [Guercil pOOOl )]. It is well understood 



that STAP techniques can improve slow-moving target detection through 
better mainlobe clutter suppression, provide better detection in combined 
clutter and jamming environments, and offer a significant increase in out- 
put signal to- interference-plus- noise-ratio (SINR). Moreover, it is also well 
understood that clutter and jamming signals often reside in a low-rank sig- 
nal subspace, which is typically much lower than the number of degrees of 
freedom of the array and the associated space-time beamformer. Due to 
the large computational complexity of the matrix inversion operation, the 
optimum STAP processor is prohibitive for practical implementation. In 
addition, another very challenging issue that is encountered by the optimal 
STAP technique is when the number of elements M in the spatio-temporal 
beamformer is large. It is well known that K > 2M independent and iden- 
tically distributed (i.i.d) training sampl es are required for the beamformer 
to achieve the steady-state performance HavkinI (|2002l )]. Thus, in dynamic 
scenarios the optimal STAP with large M usually fails or provides poor 
performance in tracking target signals contaminated by interference and 
noise. 



1 
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In the recent years, a number of innovative space-time beam- 
forming algorithms have been reported in the hterature for clutter 
and interference mitigation in radar systems. These algorithms i n- 



clude low-rank and reduced-d i mensi on techniques Haimovich I (jl991r )] 



de Lamare and Sampaio-Neto I (|2009f )]. which employ a two-stage process 



ing framework to exploit the low-rank property of the clutter and the jam- 
ming signals. The first stage performs dimensionality reduction and is 
followed by a second stage that employs a beamforming algorithm with a 
reduced dimensional filter. Another class of important space-time beam- 
forming algorithms adopt the strategy of compressive sensing and sparsity- 
aware algorithms, which exploit the fact that space-time beamformers do 
not need all their degrees of freedom to mitigate clutter and jamming sig- 
nals. These algorithms compute sparse space-time beamformers which can 
converge faster and are effective for STAP in radar systems. By exploiting 
the low-rank properties of the interference and devising sparse STAP algo- 
rithms, designers make use of prior knowledge about the clutter and the 
jamming signals. It has been recently shown that it is beneficial in terms of 
performance to also exploit prior knowledge about the environment and the 
data in the form of a known covariance data matrix. The class of space-time 
beamforming algorithms that exploit different forms of prior knowledge are 
called knowledge-aided STAP (KA-STAP) algorithms. 

The goal of this chapter is to review the recent work and advances 
in the area of space-time beamforming algorithms and t heir app l icatio n 



to radar systems. These systems include phased-arra v iMelvin I (l2004r)1 



and m ulti-input multi-output (MIMO) radar systems Haimovich et 



(l2008f)1. r nono- static and bi-static radar systems and other configurations 



Melvin I (|2004f )]. Furthermore, this chapter also describes in detail some of 



the most successful space-time beamforming algorithms that exploit low- 
rank and sparsity properties as well as the use of prior-knowledge to improve 
the performance of STAP algorithms in radar systems. 

The chapter is structured as follows. Section 11.21 describes the radar 
system under consideration and the signal model used to mathematically 
describe it. Section 11.31 formulates the problem of designing space-time 
beamformers and reviews conventional space-time beamforming algorithms. 
Section Fl .4l examines low-rank space-time beamforming algorithms, whereas 
Section 11.51 explores the concept of sparsity-aware space-time beamform- 
ing algorithms. Section 11.61 studies knowledge- aided beamforming algo- 
rithms and discusses how these techniques can be adopted in existing radar 
systems. Section 11.71 is devoted to the presentation of simulation results. 
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discussions and the comparison of a number of existing algorithms. The 
chapter ends with Section 11.81 which gives the concluding remarks of this 
chapter. 



1.2 System and Signal Models 



The system under consideration is a pulsed Doppler radar residing on an 
airborne platform. The radar antenna is a uniformly spaced linear antenna 
array consisting of N elements. The radar returns are collected in a coherent 
processing interval (CPI), which is referred to as the 3-D radar datacube 
shown in Fig. 11.11 (a), where K denotes the number of samples collected 
to cover the range interval. The data is then processed at one range of 
interest, which corresponds to a slice of the CPI datacube. This slice is 
a J X iV matrix which consists of x 1 spatial snapshots for J pulses at 
the range of interest. It is convenient to stack the matrix column- wise to 
form the M x 1 vector r(i), ter med the i-th ran ge gate spacetime snapshot, 
where M = JN aj\d l< i < K iKlemml (Eooi)]. 
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Fig. 1.1 (a) The Radar CPI datacube. (b) The STAP schematic. 



The objective of a radar is to ascertain whether targets are present in 
the data. Thus, given a space-time snapshot, radar detection is a binary 
hypothesis problem, where hypothesis corresponds to the absence of 
a target and hypothesis Hi corresponds to the presence of a target. The 
radar space-time snapshot is then expressed for each of the two hypotheses 
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in the following form 

Hn : r{i) = v{i); 

W V 7, ^^^^ 
HI : r{i) — as + v{i); 

where a is a zero-mean complex Gaussian random variable with variance 
(Tg , v{i) = rc{i) + 'r'j{i) + n{i) contains the input interference-plus-noise 
vector which consists of the clutter rdi), the jamming signal rj{i) and the 
complex white Gaussian noise These three components are assumed 
to be mutually uncorrelated. Thus, the M x M covariance matrix R of the 
undesired clutter-plus-jammer-plus- noise component can be modelled as 

R = Rc -\- Rj + Rn (l-^) 

where (•)^ represents the Hermitian transpose and E[-] denotes expecta- 
tion. The noise covariance matrix is given by i?„ = E[n{i)n^ (i)] = a^-^I , 
where is the variance of the noise and I is an identity matrix. The 
clutter signal can be modeled as the superposition of a large number of 
independent clutter patches which are evenly distributed in azimuth about 
the receiver. Thus, the clutter covariance matrix can be expressed as 

Re = E[r,r^] = J2I1 Ck,iM^h)Hr,'/) ® [a{wlMK"), (1-3) 

fe=i 1=1 

where Nr denotes the number of range ambiguities and A^c denotes the num- 
ber of clutter patches. The quantity ^ is the power of the reflected signal 
by the k, l-th clutter patch. The symbol <S> denotes Kronecker product, and 
the quantities ; and a(t:^^ ;) denote the spatial steering vector with the 
spatial frequency ; and the temporal steering vector with the normalized 
Doppler frequency ^ for the k, l-th clutter patch, respectively, which can 
be expressed as follows 



m,i) = 



g-j27r2i? 



-j27r(JV-l)?? 



1 

g-j27rro 
g-j27r2ro 



-j27r(JV-l)ro 



(1.4) 



where = jCos{(j))sin{9) and t? = fd/fr, A is the wavelength, d is the 
inter-element spacing which is normally set to half wavelength, and (p and 
6 are the elevation and the azimuth angles, respectively. The quantities 
fd and fr are the Doppler frequency and the pulse repetition frequency. 



February 12, 2013 2:27 World Scientific Bool< - 9in x 6in 



Adaptive Space-Time Beamf arming in Radar Systems 5 

respectively. The jamming covariance matrix Rj = E[rj{i)r^ [i)] can be 
written as 

= E ® Ik, (1.5) 

where ^ is the power of the q-th jammer. The vector is the spatial 

steering vector with the spatial frequency i?^ of the q-th jammer and Nj 
is the number of jamming signals. The vector s is the M x 1 normalized 
space-time steering vector in the space-time look-direction, which can be 
defined as 



/^tb{^t)a{mt), (1.6) 

where a{wt) is the K x 1 normalized temporal steering vector at the target 
Doppler frequency Wt and 6(i?t) is the A'^ x 1 normalized spatial steering 
vector in the direction provided by the target spatial frequency i?i and 
denotes the power of the target. 

1.3 Conventional Beamforming Algorithms 

In order to detect the presence of targets, each range bin is processed by 
an adaptive space-time beamformer, which is typically designed to achieve 
maximum output SINR, followed by a hypothesis test to determine the tar- 
get presence or absence. The secondary data r(i) arc taken from training 
samples, which should b e idea lly i.i.d. training samples but arc often non- 



heterogeneous [Klemml (|2002|) ]. The optimum full-rank STAP that maxi- 
mizes the SINR can obtained by solving the following minimum variance 
distortionless response (MVDR) constrained optimization given by: 

Wont — aTgiJiinw^ Rw subject to s = 1, (1-7) 

w 

where the optimal space-time MVDR beamformer lOopt is designed to max- 
imize the SINR and to maintain a normalized response in the target spatial- 
Doppler look-direction. The solution to the optimization problem above is 
described by: 

The space-time beamformer lOopt can be computed by using the above so- 
lution. Alternatively , the space-tim e beamformer can be estimated using 



adaptive algorithms HavkinI (|2Q02r )]. These algorithms include the least 
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mean-square (LMS), the conjugate gradient (CG) and the recursive least- 
squares (RLS) techniques. The computational complexity of these algo- 
rithms ranges from linear with M for the LMS to quadratic with AI for the 
CG and RLS algorithms. A common problem with the conventional adap- 
tive algorithms is that the laws that govern their convergence and tracking 
behaviors imply that they depend on M and on the eigenvalue spread of 
R. This means that their performance may degrade significantly when the 
space-time beamformer has many parameters for adaptation, which makes 
the computation of the parameters of the beamformer slow and costly. This 
problem can be addressed by some recent techniques reported in the liter- 
ature, namely, low-rank, sparsity-aware and knowledged-aided algorithms. 



1.4 Low-Rank Beamforming Algorithms 



Low-rank adaptive signal processing has been considered as a key technique 
for dealing with large systems in the last decade. The basic idea of the low- 
rank algorithms is to reduce the number of adaptive coefficients by project- 
ing the received vectors onto a lower dimensional subspace which consists of 
a set of basis vectors. The adaptation of the low-order filter within the lower 
dimensional subspace results in significant computational savings, faster 
convergence speed and better tracking performance. The first statistical 
low-rank method was based on a prin cipal-components ( PC) decomposition 



of the target- free covariance matrix Haimovich I (jl99l[ )] . Another class of 



eigen-decomposition method s was based on the cross-spectral metric (CSM) 
Goldstein and Reedl |l997cl lal)]. Both the PC and the CSM algorithms re- 



quire a high computational cost due to the eigen-decomposition. A family 
of the Krylov subspace methods has been investigated thoroughly in the re- 
cent years. Thi s class of low-rank alg o rithms include s the multistage Wiene r 
filter (MS WF) (Golstein et d\ (Il998l ) ; ICuercil ^200(1) ; ICau and Reedl (Il998h ] 
which projects the observation data onto a l ower-dimensional Krylqv sub - 
space, and the auxiliary-vect o r filters (AVF ' j iPados and Batalama I (|l999l ): 
Pados and KarvstinosI (|2Q01l ): iPados et all (|2007l) ]. These methods are rel- 
atively complex to implement in practice and may suffer from numerical 
problems despite their improved convergence and tracking performance. 
The joint domain localized (JDL) approach, which is a bea mspace reduced- 



Wang and Cai 



dimen sion algorithm, was proposed by Wang and Cai 
(|1994| )] and i nvestigated i n both homogeneous and nonhomogeneous envi- 
ronments in Adve et al. ( 2000alfbh ]. respectively. Recently, reduced-rank 
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adaptive proces sing algorithms based on joint iterativ e optimization o f 



adaptive filters de 


^amare and Samt 


>aio-Neto ( 2007a k Fa et al. 


de Lamare 


( 


20081): 


de Lamare et al. 


(120101): Fa and de Lamare 



2008); 



201l|)] 



and bas ed on an adaptive diversity-comb i ned de c imatio n and interpolatio n 
scheme Ide Lamare and Samoaio-Neto I (|2007bi boool ): iFa et al. I |2010l ): 



de Lamare et al. 



( 201ll) ] were proposed, respectively. 



M X 1 

Received signal 



Dimensionality 
Reduction 

M xD 



Dxl 



Low-Rank 
Processing 

Wo 

D X 1 



Fig. 1.2 Low-rank signal processing scheme. 



The basic idea of low-rank algorithms is to reduce the number of adap- 
tive coefficients by projecting the received vectors onto a lower dimensional 
subspace. Let So denote the M x D rank-reduction matrix with column 
vectors which form an M x 1 basis for a D-dimensional subspace, where 
D < M. Thus, the received signal r{i) is transformed into its reduced-rank 
version r_D(i) given by 

roii)^SZr{i) (1.9) 

The low-rank signal is processed by an adaptive low-rank space-time beam- 
former wu with D coefficients. This is illustrated in Fig. 2. Subse- 
quently, the decision is made based on the output of the beamformer 
y{i) = w^S^r{i). A designer can compute the parameters of the beam- 
former by solving the following constrained optimization problem: 

WD,opt = argminw^S^RSDWD subject to w^S^s = 1, (1.10) 

■WD 

The optimal low-rank MVDR solution for the above problem is given by 

s" Sd{SpRSd) SjjS sf,Rjj sd 

where Rd = S^RSd denotes the low-rank covariance matrix and sjj = 
S'j^s denotes the low-rank steering vector. The key challenge in the design 
of low-rank STAP algorithms is to find a cost-effective method to compute 
the rank-reduction matrix Sd- 
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1.4.1 Eigenvalue- decomposition-based algorithms 

The eigenvalue-decomposition (EVD)-based beamforming algorithms are 
also known as PC-based algorithms and have been originally reported as 
the eigencanceler method. These PC-based algorithms refer to the beam- 
formers constructed with a subset of the eigenvectors of the interference- 
only covariance matrix associated with the eigenvalues of largest magni- 
tu de. The first appli cation of this method to radar systems was reported 



Haimovich I (|l99ll )]. 



The basic idea of the EVD-based beamformer is to approximate the 
M X M covariance matrix R of the received data as follows: 



D 



(1.12) 



where the M x 1 vector Vd corresponds to the dth eigenvector of R and is 
the dth eigenvalue of R. By assuming that the eigenvalues are obtained in 
decreasing order of magnitude, the EVD-based method approximates R us- 
ing its D dominant eigenvectors. The rank-reduction matrix is constructed 
by using the D dominant eigenvectors as described by 



Sd ^[vi V2 ... Vd] 
The low-rank MVDR solution for the above problem is given by 



Wd = 



S d{S^RS d) ^S^s 
^"Sd{S'^RSd)-'S^s ^ s^(EtiA^'^rfO^' 



(1.13) 



(1.14) 



The EVD-based low-rank MVDR space-time beamformer described above 
does not take into account the target steering vector s when selecting a 
suitable subspace representation of the interference. Clearly, this low-rank 
space-time beamformer requires the compu tation of an EVD, which h as 
a computational cost that is cubic with M Golub and van Loan ( 20021 )]. 
In order to reduce this computational complexity, a designer can resort 
to sub s pace t racking a l gorith ms which bring the cost down to O(M^) 



Yang I (|l995l) ; iBadeau I (|2005[ )] . Another technique associated with EVD- 



based beamforming that can improve the performance of low-rank MVDR 
space-time beamformers is t he method called cross-spectral metric (CSM) 



Goldstein and Reedl (|l997c() ]. The CSM approach chooses the set of D 



eigenvectors for the rank-reduction matrix which optimizes the desired cri- 
terion, namely, the maximization of the SINR, in opposition to the PC 
method which always chooses the dominant eigenvectors. 
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1.4.2 Krylov sub space-has ed algorithms 

The first Krylov methods, namely, the conjugate gradie nt (CG) metho d 



Hestenes and Stiefel I (jl952l )] and the Lanczos algorithm Lanczos I (jl952l )] 



have been originally proposed for solving large systems of linear equations. 
These algorithms used in numerical linear algebra are mathematically iden- 
tical to each other and have been derived for Hermitian and positive def- 
inite system matrices. Other techniques have been repo rted for solving 
these problems and the Arnoldi algorithm Arnoldi I (jl95l[ )] is a computa- 
tionally efficient procedure for arbit rarily invertible syst e m mat rices. The 
multistage Wiener filter (MSW F) IColdstein and Reedl |l997cl)] and the 
auxiliary vector filtering (AVF) Pados and Batalama HJ999)] algorithms 
are based on a multistage decomposition of the linear MMSE estimator. A 
key feature of these methods is that they do not require an EVD and have a 
very good performance. It turns out that Krylov subspace algorithms that 
are used for solving very large and sparse systems of linear equations, are 
highly suitable alternatives for designing low-rank space-time beamforming 
algorithms in radar systems. The basic idea of Krylov subspace algorithms 
is to construct the rank-reduction matrix S^j with the following structure: 



SD=[qRq ■■■ R^'-'q] 



(1.15) 



where q 



Mi\\ 



and 



denotes the Euclidean norm (or the 2-norm) 



of a vector. In order to compute the basis vectors of the Krylov subspace 
(the vectors of So), a designer can either directly employ the expression 
in (jl.lSp or resort to mo re sophisticated approaches such as the Arnoldi 



Arnoldi I (|l95l[ )]. The low-rank MVDR solution for the space- 



iteration 

time beamformer using the Krylov subspace is given by 



Wd 



{S^RSd)-'S^s 



Sd{S^RS 



(1.16) 



An appealing feature of the Krylov subspace algorithms is that the required 
model order D does not scale with the system size. Indeed, when M goes to 
infinity the required D remains a finite and rel atively small value. This re- 
sult was established in [Xiao and Honig I (|2005l )]. Among the disadvantages 



of Krylov subspace methods are the relatively high computational cost of 
constructing So {0{DA'P)), the numerical instability of some implemen- 
tations and the lack of flexibility for imposing constraints on the design of 
the basis vectors. 
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1.4.3 JlO-based algorithms 

The aim of this part is to introduce the reader to low-rank beamforming 
algorithms based on joint iterative optimization (JIO) techniques. The idea 
of these methods is to design the main components of a low-rank space-time 
beamforming scheme via a general optimiz ation approach. The basic ideas 



of JIO t echniques have be en reported i n Ide Lamare and Sampaio-Neto . 
( 2007al): IPa et al\ J2OO8I): Ide Lamare I (|2008l ): Ide Lamare et al. I |20ld ): 



Fa and de Lamare I (|201ll )]. Amongst the advantages of JIO techniques 
are the flexibility to choose the optimisation algorithm and to impose con- 
straints, which provides a significant advantage over eigen-based and Krylov 
subspace methods. One disadvantage that is shared amongst the JIO tech- 
niques, eigen-based and Krylov subspace methods are the complexity asso- 
ciated with the design of the matrix S d ■ For instance, if we are to design a 
beamforming algorithm with a very large M, we still have the problem of 
having to design an M x D rank-reduction matrix S d . 

In the framework of JIO techniques, the design of the matrix S d and the 
beamforming vector wd for a fixed model order D will be dictated by the 
optimization problem and the algorithm chosen to compute the solution. 
To this end, we will focus on a generic S£) — [si S2 ... in which 

the basis vectors s^, d — 1,2, . . . , D will be obtained via an optimization 
algorithm and iterations between the S^i and wd will be performed. The 
JIO method consists of solving the following optimization problem 



['S'_D,opt,i^_D,opt] = arg min w%S"RSdWd, 

x(i) 



subject to w^S^s = 1 



(1.17) 



where it should be remarked that the optimization problem in (jl.l7p is non 
convex, however, the algorithms do not present convergence problems. Nu- 
merical studies with JIO methods indicate that the minima are identical and 
global. Proofs of global convergence have been established wit h different 



versio n s of JIO schemes Ide Lamare and Sampaio-Neto I (l2007al): iFa et al. 
( 2008 ): de Lamare ( 20081) : de Lamare et al. I ( 2010l ): Fa and de Lamare 



201 ll )] , which demonstrate that a least-squares (LS) algorithm converges 



to the reduced-rank Wiener filter. 

In order to solve the ab ove problem, we resort to the method of Lagrange 
multipliers HavkinI (|2002[ )] and transform the constrained optimization into 
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an unconstrained one expressed by the Lagrangian 



C{Sd,wd) = wIsIRSdwd + Kwd^dS - 1) 



(1.18) 



where A is a scalar Lagrange multipher. By fixing wo, minimizing p.l^ 
with respect to Sd and solving for A, we obtain 



R 



,"R^^wdsHR-^s' 



(1.19) 



where R — E[r{i)r^ [i)] and Rw — E[w dw^]. By fixing Sdi minimizing 
(|1.18[) with respect to wd and solving for A, we arrive at the expression 



R~ 



s^R- 



(1.20) 



where Rd = E[Sar{i)r" {i)SD] -B[r£,(i)rg(i)], sd = S"s. Note that 
the expressions in (|1.19p and (|1.20p are not closed- form solutions for wjj 
and Sd since (I1.19P is a function of wd and (|1.20l) depends on Sd- Thus, 
it is necessary to iterate (|1.19l) and (|1.20p with initi al values to obtain a 
solution . Unlike the Kr ylov subspace-based me t hods [Goldstein and Reed 
(|l997ah ] and the AVF |Pados and KarvstinosI |200ll )] methods, the JIO 
scheme provides an iterative exchange of information between the low-rank 
beamformer and the rank-reduction matrix and leads to a simpler adaptive 
implementation. The key strategy lies in the joint optimization of the filters. 
The rank D must be set by the designer to ensure appropriate performance 
or can be estimated via another algorithm. In terms of complexity, the JIO 
techniques have a computational cost that is related to the optimization 
algorithm. With recursive LS algorithms the complexity is quadratic with 
M ((0(M^)), whereas the complexity can b e as low as linear with M when 



stocha stic gradient algorithms are adopted jde Lamare and Sampaio-Neto 
|2009h ]. 



1.4.4 JIDF-based algorithms 

This section reviews a low-rank space-time beamforming technique based 
on the joint interpolation, decimation a n d filtering (JID F) concept 



de Lamare and Sampaio-Neto I |2007bll2009l ) : iFa et al. I |2010l )]. The JIDF 



approach allows a designer to compute the parameters of the rank- reduction 
matrix and the low-rank space-time beamformer with a low complexity. 
The motivation for designing a rank-reduction matrix based on interpola- 
tion and decimation comes from two observations. The first is that rank 
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reduction can be performed by constructing new samples with interpola- 
tors and eliminating (decimating) samples that are not useful in the STAP 
design. The second comes from the structure of the rank-reduction matrix, 
whose columns are a set of vectors formed by the interpolators and the 
decimators. 

In the JIDF scheme, the number of elements for adaptive processing 
is substantially reduced, resulting in considerable computational savings 
and very fast convergence performance for the radar applications. The 
M X 1 received vector r{i) is processed by a multiple processing branch 
(MPB) scheme with B branches, where each spatio-temporal processing 
branch contains an interpolator, a decimation unit and a low-rank space- 
time beamformer. In the 6-th branch, the received vector r(i) is filtered by 
the interpolator Vb = [vb,o Vb,i ■ . ■ Vb^i-i]'^ with / coefficients, resulting in 
an interpolated received vector r(,(i) with M samples, which is expressed 
by 

nil) = r(z), (1.21) 
where the M x M Toeplitz convolution matrix is given by 

... 



Vb,0 
VbJ-l 








Vb,I-l 






"6,0 



(1.22) 



The vector rb{i) can be expressed in an alternative way that is useful for 
the design of the JIDF scheme and is described by 

nil) = V^r{i) = Ua{i)vb, (1.23) 

where the M x I matrix Sto{i) with the samples of r{i) has a Hankel struc- 
ture and is described by 

ro(?:) ri{i) ...ri_i(i)' 
ri{i) r2{i) ... n(«) 



5Ro(«) 



rM-2{i) rM-i{i) 
.rM-i{i) 



(1.24) 
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(1.25) 



The dimensionality reduction is performed by a decimation unit with DxM 
decimation matrices Db that transforms r/(i) into D x 1 vectors rD,b{i) 
with b = 1, . . . ,B, where D — AI/L is the rank of the resuhing system of 
equations that will be generated and L is the decimation factor. The 13 x 1 
vector r]j^h{i) for branch b is expressed by 

rD,b - ^g,fcr(z) = DD,bV^r{i) 

where SD,b is the rank-reduction matrix and the vector r£),b(z) for branch 
b is used in the minimization of the output power for branch b. The output 
at the end of the JIDF scheme y{i) is selected according to 

yii) = VbAi) when 6^ = argmin ly^p (1-26) 

b 

where i? is a parameter to be set by the designer. For the computation of 
the parameters of the JIDF scheme, it is fundamental to express the output 
yb{i) as a function of the interpolator Vh, the decimation matrix Dd^i, and 
the low-rank space-time beamformer wjj^i, as follows: 

yb{i) = w^j,S"^br{i) 

= w",^DD,b^oii)v, 

where the expression (|1.27p indicates that the dimensionality reduction car- 
ried out by the JIDF scheme depends on finding appropriate ■»(,, DD,b and 
WD,b- Unlike the previously discussed low-rank beamforming techniques, 
the JIDF is able to substantially reduce the cost of the rank-reduction ma- 
trix. 

The parameters of the JIDF scheme that perform low-rank space-time 
MVDR beamforming can be computed by solving the following optimiza- 
tion problem 

[iOD,opt, ■"opt, -C>D,bJ = arg min w" ^E[DD,b^o{i)vv"^" {i)D'j^ ^]wd 
subject to w^j^Do.bSoVb — 1, 

(1.28) 

where So is M x I steering matrix with a Hankel structure, which has the 
same form as 5io(i) and is given by 

so(«) si{i) . . . si_i(i) 
Sl(j) S2(«) ••■ si(j) 



(1.27) 



So{^) 



Sj\/-2(j) SAf-l(j) ... 
SM-l(j) ... 



(1.29) 
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The constrained optimization in (|1.28|) can be transformed into an uncon- 
strained optimization problem by using the method of Lagrange multipliers, 
which results in 

£{wD,b,Vb,DD.b) = w'^E[DDMi)vv"U^ {i)D'^ .jwD+Kw'^DD^bSoV-l), 

(1.30) 

where A is a Lagrange multiplier. 

The strategy to compute the parameters of the low-rank space-time 
beamformer based on the JIDF scheme is to minimize the cost function 
with respect to a set of parameters and fix the remaining parameters. By 
minimizing (|1.30p with respect to v^, we obtain 

R-l„ 

v,b''v,b 
''v,b-^v,b''v,b 

where Rv,b — £'[r^,br^j] is the I x I autocorrelation matrix, r^^j, = 
-D^ ^R^WBfi, and s^.b — i^S^wdj, is the Ix low-rank steering vector. 
By minimizing (jl.30p with respect to wo.b^ we have 

Rw,b^it',b 

WD,b= 1 , (1-32) 

^w,b^w,b^w,b 

where Rw,b — Elr^^fwi^ '^^ D x D autocorrelation matrix, r^^b = 
DD,bRoVb, and Sw^b — Do.bSoVb is the D x I low-rank steering vector. In 
order to compute Vb and Wjj bi a designer needs to iterate them for each 
processing branch b. 

The decimation matrix b is selected to minimize the square of the 
output of the beamformer yb{i) obtained for all the B branches 

Do.b^Do.bM when 6s arg min \yb{i)\'^ , (1.33) 

l<b<B 

The design of the decimation matrix D u b imposes constraints on the values 
of the elements of the matrix such that they only take the value zero or 
one. Since the optimal approach for the design of Djj b corresponds to an 
exhaustive search, we consider a suboptimal technique that employs pre- 
stored patterns. The decimation scheme employs a structure formed in the 
following way 

So.b = [</'64 0b,2 Am], (1-34) 
where <pb d an M x 1 vector composed of a single one and zeros elsewhere 
as described by 

ct>b^^ = 1, 0^_^], (1.35) 
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where Zt^d is the number of zeros before the only element equal to one. We 
set the value of zt^d in a deterministic way which can be expressed as 

zt,d = ^^(d-l) + ib-l). (1.36) 

It is necessary to iterate (|1.3ip . (I1.32[) and (ll.33|) in an alternated form 
(one followed by the other) with an initial value to obtain a solution. The 
expectations can be estimated either via time averages or by instantaneous 
estimates and with the help of adaptive algorithms. 



1.5 Sparsity- Aware Beamforming Algorithms 

This section considers space-time beamforming algorithms that exploit the 
sparsity encountered in the data processed by radar systems. In particu- 
lar, the motivation for exploiting the sparsity of data vectors observed by 
radar systems is given and a brief discussion on the suitability of sparsity- 
aware algorithms for radar applications is provided. A general approach 
to design space-time beamforming algorithms based on the Zi-norm regu- 
larization is described. The main principle is to employ a reduced number 
of weights to suppress the clutter and the jamming signals encountered in 
radar applications. 

Recently, motivated by compressive sensing (CS) techniques used in 
radar, several authors have consi dered CS ideas for mov i ng target indicatio n 
(MTI) and STAP problems [Mar ia and Fuchs I |2006l )]- |Selesnick I |2O10l )]. 



The core notion in CS is to regularize a linear inverse p roblem by including 



prior knowledge that the signal of interest is sparse Parker and Potter 



(|2010f )]. These works on space-time beamforming techniques based on CS 
rely on the recovery of the clutter power in angle-Doppler plane, which is 
usually carried out via convex optimization tools. However, these methods 
are based on linear programming and have a quite high computational 
complexity {0{K^)), where K is the dimension of the angle-Dopplcr plane. 
In this section, we describe the concept of a sparsity-aware STAP (SA- 
STAP) algorithm that can improve the detection capability using a small 
number of snapshot. To overcome the high complexity of the CS-STAP 
type algorithm, we design the STAP algorithm with another strategy, by 
imposing the sparse regularization to the minimum variance (MV) cost 
function. Since the interference variance has often a low-rank property, we 
assume that a number of samples of the data cube are not meaningful for 
processing and the optimal STAP beamformer is sparse, or nearly sparse. 
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Then, we exploit this feature by using a Zi-norm regularization. With this 
motivation, the STAP algorithm design becomes a mixed /i-norm and I2- 
norm optimization problem. 

The conventional space-time beamforming algorithms do not exploit 
the sparsity of the received signals. In this exposition, it is assumed that 
a number of samples of the data cube are not meaningful for processing 
and a reduced number of active weights of the space-time beamformer can 
effectively suppress the clutter and the jamming signals. Specifically, a 
sparse regularization is imposed to the space-time MVDR beamforming 
design. Thus, the space-time beamformer design can be described as the 
following optimization problem 

Wont — avgimnw^Rw 

- (1.37) 
subject to w s = 1 and \ \w\\i = 0, 

where the objective of the Zi-norm regula rization is to force th e comp onents 



of the space-time beamformer w to zero jAngelosante et al. I pOlOf l]. This 
problem can be solved using the method of Lagrange multipliers, which 
results in the following unconstrained cost function 

C{w,a,X) = w" Rw + a{w" s - 1) + X{\\w\\i), (1.38) 

The unconstrained cost function above is convex, however, it is non- 
differentiable which makes it difficult for one to use the method of La- 
grange Multipliers directly and obtain an expression for the space-time 
beamformer. To this end, the following approximation to the regulariza- 
tion term is employed 

\\w\\i^w"Aw, (1.39) 

where 

A - diag f . . . I \—) , (1.40) 

where e is a small positive constant. Simultaneously, we assume that the 
partial derivative of Aw with respect to w* is given by 

-g^^Au,. (1.41) 

With the development above, an approximation to the unconstrained cost 
function can be employed as described by 

C{w, a, A) « w"Rw + a{w"s - 1) + Xw" Aw, (1.42) 
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By computing the gradient terms with respect to w* and a and equat- 
ing them to zero, we obtain the following expression for the space-time 
beamformer 

w = ^— /. , , . 1.43 

Comparing (jl.43l) with the conventional optimal space-time beamformer 
in (|1.8|) . we find that there is an additional term A A in the inverse of the 
interference covariance matrix ii, which is due to the Zi-norm regulariza- 
tion. The term A is a positive scalar which provides a trade-off between 
the sparsity and the output interference power. The larger the chosen A, 
the more components are shrunk to zero [Zibulevskv and Elad ] lioio)]. It 
should also be remarked that the expression for the beamformer in ()1.43p 
is not a closed-form solution since A is a function of w. Thus it is neces- 
sary to develop an iterative procedure to compute the parameters of the 
space-time beamformer. 



1.6 Knowledge- Aided Beamforming Algorithms 



Although STAP techniques are considered efficient tools for detection of 
slow targets b y airborne radar systems in strong clutter environments 
Klemml(|2003 l]. due to the very large number of degrees of freedom (DoFs) 



conventional space-time beamformers have a slow convergence and re- 
quire about twice the DoFs of the independent and identically distributed 
(IIP) training snap shots to yield an average performance loss of roughly 
3dB Ward I ([l99J)]. In real scenarios, it is hard to obtain so many IID 
training snaps h ots, especially in heterogeneou s envir onments. Low-rank 
Guerci ( "2000[l1- de Lamare and Sampaio-Neto ( 20091 )] and sparsity-aware 



Maria and Fuchs 



( 20061 1]- [Yang et al. I (|201lh ] methods have been consid- 



ered to counteract the slow convergence of the conventional space-time 
beamformers. Nevertheless, there are other alternatives to improve the 
training of STAP algorithms and improve their performance. These other 
methods can also be combined with the techniques previously discussed. 
Recently developed knowledge-aided (KA) STAP algorithms have received 
a growing interest and become a key conce p t for the next ge neration of 



a growing interest ana pecome a Jtey conce p t tor tne next ge neration ot 
adaptive radar systems Wicks et al. I ( 20061 )]- Fa et al. ( 2O10l )]. The core 
idea of KA-STAP is to incorporate prior knowledge, provided by digital 
elevation maps, land cover databases, road maps, the Global Position- 
ing System (GPS), previous scanning data and other known features, to 
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compute estimates of the clutter covariance matrix with high accuracy 



Melvin and Guerci I (j2006t )]. Prior work on KA-STAP algorithms include 



the exploitati on of prior knowledge of the clutter ridge to form the STAP 
filter weights Melvin and Showman ( 2006I )]. use of prior knowledge about 



the terrain [Capraro et al. I (|2006[ )] and prior knowle dge about th e covar i 



ance matr i x of t he clutter and the jamming signals [Blunt et al. ( 20061) ]- 



Fa et al. 



In this section, we discuss a strategy to mitigate the deleterious effects 
of the heterogeneity in the secondary data, which makes use of a priori 
knowledge of the clutter covari ance matrix and has r ecently gained si gnifi- 



cant attention in the literature (Wicks et al7\ (|2006[ )]- |Fa et al. I |2010l )]. In 



KA-STAP techniques, there are two basic tasks that need to be addressed: 
the first one is how to obtain prior knowledge from the terrain knowledge 
of the clutter and how to estimate the real interference covariance matrix 



with the prior knowledge [Wicks et al. I (|2006l )]- [Capraro et al. I (|2006l )] and 
the second is how to apply the covariance matrix estimates in the design 



of the space-time beamforming algorithm Blunt et al. I (|2006r )]- Fa et al. 
(|2010l )]. We first review how a designer can obtain prior knowledge of the 



clutter and employ this knowledge to build a known covariance matrix i?o ■ 
Then, we present a method to combine this prior knowledge with com- 
monly used estimation techniques to compute the covariance matrix of the 
received vector r{i), resulting in a combined covariance matrix estimate 
Rc for use in the space-time beamformer that is more accurate and has an 
enhanced performance. 

The optimal space-time beamformer employs the following expression 
to compute its parameters 



R 



s"R 's' 



(1.44) 



where an estimate of the covariance matrix is typically obtained by 



(1.45) 



k=l 



where r{k) is taken from secondary data. The e stimate R can be sufficientl y 
accurate when K is at least twice as great as M Brennan and Reed I (|l973r )] 
and the training samples are assumed i.i.d. However, it is by now well 
understood that the clutter environments are often heterogeneous and this 
leads to performance degradation on space-time beamforming. KA-STAP 
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techni ques can significantly help to combat the heterogeneity [Stoica et al. 



torn ]. 

With KA techniques the covariance matrix is estimated by combin- 
ing an initial guess of the covariance matrix Ro derived from the digital 
terrain database or the data probed by radar in previous scans, and the 
sample average covariance matrix estimate in the present scan R so that 

Rc = aRo + (1 - a)R, (1.46) 

where < a < 1 . Alternatively, this principle can be applied to the inverse 
of the covariance matrix estimate 

R-' ^ fjR-' + {1 - rj)R~\ (1.47) 

where < < 1. 

In order to compute the parameter 77, we need to consider the optimiza- 
tion problem 

Vopt = argmin Rw, (1-48) 

where we use the relation 

w{i) = 1TW0 + {1 - v)w{i), (1.49) 
where w{i) — R^s, Wo — Rq^^ and w — R s. We can obtain the 



optimal value for rj by ec 
to zero, which results in 



uating the gr a dient of the cost function in (jl.48p 



Stoica et al. I (|2008[ )] 



^{s"{R ^ R-^)RR ^s) 



Vopt = ^ / _i ° \ (1-50) 

s"{R-^-R )R{R-^ R )s 

Since R above is unknown, we have to estimate it in real time using either 
time averages or adaptive algorithms. 



1.7 Simulations 

In this section, we assess the performance of the space-time beamforming 
algorithms discussed in this chapter using simulated radar data. Specifi- 
cally, we consider the optimal MVDR space-time beamforming algorithm 
that assumes perfect knowledge of the covariance matrix of the received 
data, and the MVDR space-time beamformer using the sample matrix in- 
version (SMI-MVDR). The low-rank space-time beamforming algorithms 
using EVD (LR-EVD), the Krylov subspace approach (LR-Krylov), the 
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JIO (LR-JIO) and the JIDF (LR-JIDF) algorithms are also considered with 
a rank equal to D. We also consider the sparsity-aware (SA-MVDR) and 
the knowledged-aided (KA-MVDR) space-time beamforming algorithms. 
All the analyzed algorithms estimate the statistical quantities via time- 
averages in a similar way to a least-squares method. The parameters of the 
simulated radar platform are shown in Table 11.11 For all simulations, we 
assume the presence of a mixture of two broadband jammers at —45° and 
60° with jammer-to- noise ratio (JNR) equal to 40 dB. The clutter-to-noise- 
ratio (CNR) is fixed at 40 dB. All the results presented are averages over 
1000 independent Monte-Carlo runs. 



Parameter Value 

Antenna array Sideway-looking array (SLA) 
Carrier frequency (/c) 450 MHz 

Transmit pattern Uniform 



In the first experiment, we assess the SINR performance of the different 
space-time beamforming algorithms as shown in Fig. 11.31 The algorithms 
are simulated over 800 snapshots and the SNR is set to 10 dB. The re- 
sults show that the LR-JIDF algorithm achieves the best results, followed 
by the LR-JIO, the KA-MVDR, the SA-MVDR, the LR-Krylov, the LR- 
EIG and the SMI MVDR algorithms. The curves indicate that the use of 
low-rank algorithms is highly beneficial to the performance of space-time 
beamforming algorithms in radar systems. In particular, the LR-JIDF and 
LR-JIO algorithms have a very fast convergence performance. It should 
also be remarked that the SA-MVDR and KA-MVDR algorithms obtain a 
performance that is significantly better than the conventional SMI MVDR 
algorithm. Since the SA-MVDR and KA-MVDR techniques are modifica- 
tions of the SMI MVDR techniques exploiting sparsity and prior knowledge 
about the covariance matrix, respectively, it is interesting to note that by 
exploiting these properties it is possible to significantly improve the perfor- 
mance of beamforming algorithms. 




75 m/s 
9000 m 
40 dB 



300 Hz 



8 
8 
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30 




Snapshots 

Fig. 1.3 SINR Performance against the snapshots. Parameters; B = 8, I = 8 and 
D = G, M = 64, R[0] = 0.017, e = 0.1. 



We evaluate the SINR performance against the target Doppler frequency 
at the main bean look angle for our proposed algorithms and other exist- 
ing algorithms, which are illustrated in Fig. 11.41 The potential Doppler 
frequency space form —100 to 100 Hz is examined and 100 snapshots are 
used to train the beamformers. The plots show that the analyzed algo- 
rithms converge and approach the optimum in a short time, and form a 
deep null to cancel the main beam clutter. Again, the LR-JIDF algorithm 
outperforms the other analyzed algorithms. 

In the third example, the probability of detection Pd versus SNR per- 
formance is presented for all schemes using 200 snapshots as the training 
data as shown in Fig. 11.51 The false alarm rate PFA is set to 10~^ and we 
suppose the target injected in the boresight (0°) with Doppler frequency 
lOOHz. The figure illustrates that the analyzed algorithms provide sub- 
optimal detection performance using short support data. Note that for 
Pd — 0.9 (90% percent), the LR-JIDF and LR-JIO schemes are within less 
than 1 dB from the performance of the optimal MVDR algorithm. The 
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30 




Target Doppler Frequency ^ (Hz) 



Fig. 1.4 SINR Performance against the target Doppler frequency. Parameters: B = S, 
/ = 8 and D = 6, M = 64, R[0] = 0.01/, e = 0.1. 



remaining techniques exhibit increasing performance losses as compared to 
the optimal MVDR algorithm and it should be noted that the conventional 
SMI MVDR method has a performance degradation of up to 5 dB for the 
same performance measured in terms of Pp. This suggest that the ap- 
plication of more sophistical space-time beamforming algorithms is key to 
achieving an improved performance. 

At last, we detail the computational complexity in terms of multipli- 
cations of the analyzed schemes in Fig. 11.61 The curves show that the 
computational complexity of the LR-JIDF and LR-Krylov algorithms is 
significantly lower than the remaining algorithms. Indeed, there is a sig- 
nificant computational advantage obtained by using the LR-JIDF and LR- 
Krylov algorithms and this advantage becomes more pronounced as M is 
increased. The other analyzed algorithms have a complexity that scales as 
a cubic function of M. This high complexity can be mitigated by the use of 
adaptive algorithms, which can reduce the computational cost by at least 
an order of magnitude. 
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SNR (dB) 



Fig. 1.5 Probability of detection against SNR. Parameters: B 
M = 64, R\0] = O.Oli", e = 0.1. 



and D = 6, 



1.8 Concluding Remarks 



This chapter has presented recent advances on space-time beamforming 
algorithms for phased-array radar systems and investigated their perfor- 
mance via computer simulations. Specifically, we have considered in detail 
some of the most successful space-time beamforming algorithms that ex- 
ploit low-rank and sparsity properties as well as the use of prior-knowledge 
to improve the performance of STAP algorithms. The results of our studies 
suggest that the low-rank algorithms have a substantial performance advan- 
tage over conventional MVDR space-time beamforming algorithms. Fur- 
thermore, the use of sparsity- aware and knowledge-aided strategies is also 
able to improve the performance of space-beamforming algorithms, and can 
be combined with low-rank schemes. These beamforming algorithms can 
be also applied to MIMO radar systems, mono-static and bi-static radar 
systems and other sensing applications such as sonar systems. 
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Fig. 1.6 Computational complexity in terms of multiplications of the analyzed space- 
time beamforming algorithms. Parameters: S = 8, / = 8 and _D = 6. 
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Chapter 1 

Adaptive Space-Time Beamforming in 
Radar Systems 



1.1 Introduction 



Space-time 

adaptive processing (STAP) techniques [Kleminl |2002h ]. |Melvin I |2004[ )] 



have been thoroughly investigated in the last decades as a key enabling 
technology for advanced airb orne radar applications fol lowing the seminal 



work by Brennan and Reed Brennan and Reed I (jl973h ]. A great deal of 



attention has been given to STAP algorithms and different strategies to 
design space- time beamfo r mers to mitigat e the e ffect of clutter and jam- 



ming signals [Reed et all (|1974t )]- [Guercil pOOOl )]. It is well understood 



that STAP techniques can improve slow-moving target detection through 
better mainlobe clutter suppression, provide better detection in combined 
clutter and jamming environments, and offer a significant increase in out- 
put signal to- interference-plus- noise-ratio (SINR). Moreover, it is also well 
understood that clutter and jamming signals often reside in a low-rank sig- 
nal subspace, which is typically much lower than the number of degrees of 
freedom of the array and the associated space-time beamformer. Due to 
the large computational complexity of the matrix inversion operation, the 
optimum STAP processor is prohibitive for practical implementation. In 
addition, another very challenging issue that is encountered by the optimal 
STAP technique is when the number of elements M in the spatio-temporal 
beamformer is large. It is well known that K > 2M independent and iden- 
tically distributed (i.i.d) training sampl es are required for the beamformer 
to achieve the steady-state performance HavkinI (|2002l )]. Thus, in dynamic 
scenarios the optimal STAP with large M usually fails or provides poor 
performance in tracking target signals contaminated by interference and 
noise. 
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In the recent years, a number of innovative space-time beam- 
forming algorithms have been reported in the hterature for clutter 
and interference mitigation in radar systems. These algorithms i n- 



clude low-rank and reduced-d i mensi on techniques Haimovich I (jl991r )] 



de Lamare and Sampaio-Neto I (|2009f )]. which employ a two-stage process 



ing framework to exploit the low-rank property of the clutter and the jam- 
ming signals. The first stage performs dimensionality reduction and is 
followed by a second stage that employs a beamforming algorithm with a 
reduced dimensional filter. Another class of important space-time beam- 
forming algorithms adopt the strategy of compressive sensing and sparsity- 
aware algorithms, which exploit the fact that space-time beamformers do 
not need all their degrees of freedom to mitigate clutter and jamming sig- 
nals. These algorithms compute sparse space-time beamformers which can 
converge faster and are effective for STAP in radar systems. By exploiting 
the low-rank properties of the interference and devising sparse STAP algo- 
rithms, designers make use of prior knowledge about the clutter and the 
jamming signals. It has been recently shown that it is beneficial in terms of 
performance to also exploit prior knowledge about the environment and the 
data in the form of a known covariance data matrix. The class of space-time 
beamforming algorithms that exploit different forms of prior knowledge are 
called knowledge-aided STAP (KA-STAP) algorithms. 

The goal of this chapter is to review the recent work and advances 
in the area of space-time beamforming algorithms and t heir app l icatio n 



to radar systems. These systems include phased-arra v iMelvin I (l2004r)1 



and m ulti-input multi-output (MIMO) radar systems Haimovich et 



(l2008f)1. r nono- static and bi-static radar systems and other configurations 



Melvin I (|2004f )]. Furthermore, this chapter also describes in detail some of 



the most successful space-time beamforming algorithms that exploit low- 
rank and sparsity properties as well as the use of prior-knowledge to improve 
the performance of STAP algorithms in radar systems. 

The chapter is structured as follows. Section 11.21 describes the radar 
system under consideration and the signal model used to mathematically 
describe it. Section 11.31 formulates the problem of designing space-time 
beamformers and reviews conventional space-time beamforming algorithms. 
Section Fl .4l examines low-rank space-time beamforming algorithms, whereas 
Section 11.51 explores the concept of sparsity-aware space-time beamform- 
ing algorithms. Section 11.61 studies knowledge- aided beamforming algo- 
rithms and discusses how these techniques can be adopted in existing radar 
systems. Section 11.71 is devoted to the presentation of simulation results. 
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discussions and the comparison of a number of existing algorithms. The 
chapter ends with Section 11.81 which gives the concluding remarks of this 
chapter. 



1.2 System and Signal Models 



The system under consideration is a pulsed Doppler radar residing on an 
airborne platform. The radar antenna is a uniformly spaced linear antenna 
array consisting of N elements. The radar returns are collected in a coherent 
processing interval (CPI), which is referred to as the 3-D radar datacube 
shown in Fig. 11.11 (a), where K denotes the number of samples collected 
to cover the range interval. The data is then processed at one range of 
interest, which corresponds to a slice of the CPI datacube. This slice is 
a J X iV matrix which consists of x 1 spatial snapshots for J pulses at 
the range of interest. It is convenient to stack the matrix column- wise to 
form the M x 1 vector r(i), ter med the i-th ran ge gate spacetime snapshot, 
where M = JN aj\d l< i < K iKlemml (Eooi)]. 



/-th Range 
Gates 
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Ho: Target Absent 
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(b) 



Fig. 1.1 (a) The Radar CPI datacube. (b) The STAP schematic. 



The objective of a radar is to ascertain whether targets are present in 
the data. Thus, given a space-time snapshot, radar detection is a binary 
hypothesis problem, where hypothesis corresponds to the absence of 
a target and hypothesis Hi corresponds to the presence of a target. The 
radar space-time snapshot is then expressed for each of the two hypotheses 
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in the following form 

Hn : r{i) = v{i); 

W V 7, ^^^^ 
HI : r{i) — as + v{i); 

where a is a zero-mean complex Gaussian random variable with variance 
(Tg , v{i) = rc{i) + 'r'j{i) + n{i) contains the input interference-plus-noise 
vector which consists of the clutter rdi), the jamming signal rj{i) and the 
complex white Gaussian noise These three components are assumed 
to be mutually uncorrelated. Thus, the M x M covariance matrix R of the 
undesired clutter-plus-jammer-plus- noise component can be modelled as 

R = Rc -\- Rj + Rn (l-^) 

where (•)^ represents the Hermitian transpose and E[-] denotes expecta- 
tion. The noise covariance matrix is given by i?„ = E[n{i)n^ (i)] = a^-^I , 
where is the variance of the noise and I is an identity matrix. The 
clutter signal can be modeled as the superposition of a large number of 
independent clutter patches which are evenly distributed in azimuth about 
the receiver. Thus, the clutter covariance matrix can be expressed as 

Re = E[r,r^] = J2I1 Ck,iM^h)Hr,'/) ® [a{wlMK"), (1-3) 

fe=i 1=1 

where Nr denotes the number of range ambiguities and A^c denotes the num- 
ber of clutter patches. The quantity ^ is the power of the reflected signal 
by the k, l-th clutter patch. The symbol <S> denotes Kronecker product, and 
the quantities ; and a(t:^^ ;) denote the spatial steering vector with the 
spatial frequency ; and the temporal steering vector with the normalized 
Doppler frequency ^ for the k, l-th clutter patch, respectively, which can 
be expressed as follows 



m,i) = 



g-j27r2i? 



-j27r(JV-l)?? 



1 

g-j27rro 
g-j27r2ro 



-j27r(JV-l)ro 



(1.4) 



where = jCos{(j))sin{9) and t? = fd/fr, A is the wavelength, d is the 
inter-element spacing which is normally set to half wavelength, and (p and 
6 are the elevation and the azimuth angles, respectively. The quantities 
fd and fr are the Doppler frequency and the pulse repetition frequency. 
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respectively. The jamming covariance matrix Rj = E[rj{i)r^ [i)] can be 
written as 

= E ® Ik, (1.5) 

where ^ is the power of the q-th jammer. The vector is the spatial 

steering vector with the spatial frequency i?^ of the q-th jammer and Nj 
is the number of jamming signals. The vector s is the M x 1 normalized 
space-time steering vector in the space-time look-direction, which can be 
defined as 



/^tb{^t)a{mt), (1.6) 

where a{wt) is the K x 1 normalized temporal steering vector at the target 
Doppler frequency Wt and 6(i?t) is the A'^ x 1 normalized spatial steering 
vector in the direction provided by the target spatial frequency i?i and 
denotes the power of the target. 

1.3 Conventional Beamforming Algorithms 

In order to detect the presence of targets, each range bin is processed by 
an adaptive space-time beamformer, which is typically designed to achieve 
maximum output SINR, followed by a hypothesis test to determine the tar- 
get presence or absence. The secondary data r(i) arc taken from training 
samples, which should b e idea lly i.i.d. training samples but arc often non- 



heterogeneous [Klemml (|2002|) ]. The optimum full-rank STAP that maxi- 
mizes the SINR can obtained by solving the following minimum variance 
distortionless response (MVDR) constrained optimization given by: 

Wont — aTgiJiinw^ Rw subject to s = 1, (1-7) 

w 

where the optimal space-time MVDR beamformer lOopt is designed to max- 
imize the SINR and to maintain a normalized response in the target spatial- 
Doppler look-direction. The solution to the optimization problem above is 
described by: 

The space-time beamformer lOopt can be computed by using the above so- 
lution. Alternatively , the space-tim e beamformer can be estimated using 



adaptive algorithms HavkinI (|2Q02r )]. These algorithms include the least 
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mean-square (LMS), the conjugate gradient (CG) and the recursive least- 
squares (RLS) techniques. The computational complexity of these algo- 
rithms ranges from linear with M for the LMS to quadratic with AI for the 
CG and RLS algorithms. A common problem with the conventional adap- 
tive algorithms is that the laws that govern their convergence and tracking 
behaviors imply that they depend on M and on the eigenvalue spread of 
R. This means that their performance may degrade significantly when the 
space-time beamformer has many parameters for adaptation, which makes 
the computation of the parameters of the beamformer slow and costly. This 
problem can be addressed by some recent techniques reported in the liter- 
ature, namely, low-rank, sparsity-aware and knowledged-aided algorithms. 



1.4 Low-Rank Beamforming Algorithms 



Low-rank adaptive signal processing has been considered as a key technique 
for dealing with large systems in the last decade. The basic idea of the low- 
rank algorithms is to reduce the number of adaptive coefficients by project- 
ing the received vectors onto a lower dimensional subspace which consists of 
a set of basis vectors. The adaptation of the low-order filter within the lower 
dimensional subspace results in significant computational savings, faster 
convergence speed and better tracking performance. The first statistical 
low-rank method was based on a prin cipal-components ( PC) decomposition 



of the target- free covariance matrix Haimovich I (jl99l[ )] . Another class of 



eigen-decomposition method s was based on the cross-spectral metric (CSM) 
Goldstein and Reedl |l997cl lal)]. Both the PC and the CSM algorithms re- 



quire a high computational cost due to the eigen-decomposition. A family 
of the Krylov subspace methods has been investigated thoroughly in the re- 
cent years. Thi s class of low-rank alg o rithms include s the multistage Wiene r 
filter (MS WF) (Golstein et d\ (Il998l ) ; ICuercil ^200(1) ; ICau and Reedl (Il998h ] 
which projects the observation data onto a l ower-dimensional Krylqv sub - 
space, and the auxiliary-vect o r filters (AVF ' j iPados and Batalama I (|l999l ): 
Pados and KarvstinosI (|2Q01l ): iPados et all (|2007l) ]. These methods are rel- 
atively complex to implement in practice and may suffer from numerical 
problems despite their improved convergence and tracking performance. 
The joint domain localized (JDL) approach, which is a bea mspace reduced- 



Wang and Cai 



dimen sion algorithm, was proposed by Wang and Cai 
(|1994| )] and i nvestigated i n both homogeneous and nonhomogeneous envi- 
ronments in Adve et al. ( 2000alfbh ]. respectively. Recently, reduced-rank 
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adaptive proces sing algorithms based on joint iterativ e optimization o f 



adaptive filters de 


^amare and Samt 


>aio-Neto ( 2007a k Fa et al. 


de Lamare 


( 


20081): 


de Lamare et al. 


(120101): Fa and de Lamare 



2008); 



201l|)] 



and bas ed on an adaptive diversity-comb i ned de c imatio n and interpolatio n 
scheme Ide Lamare and Samoaio-Neto I (|2007bi boool ): iFa et al. I |2010l ): 



de Lamare et al. 



( 201ll) ] were proposed, respectively. 



M X 1 

Received signal 



Dimensionality 
Reduction 

M xD 



Dxl 



Low-Rank 
Processing 

Wo 

D X 1 



Fig. 1.2 Low-rank signal processing scheme. 



The basic idea of low-rank algorithms is to reduce the number of adap- 
tive coefficients by projecting the received vectors onto a lower dimensional 
subspace. Let So denote the M x D rank-reduction matrix with column 
vectors which form an M x 1 basis for a D-dimensional subspace, where 
D < M. Thus, the received signal r{i) is transformed into its reduced-rank 
version r_D(i) given by 

roii)^SZr{i) (1.9) 

The low-rank signal is processed by an adaptive low-rank space-time beam- 
former wu with D coefficients. This is illustrated in Fig. 2. Subse- 
quently, the decision is made based on the output of the beamformer 
y{i) = w^S^r{i). A designer can compute the parameters of the beam- 
former by solving the following constrained optimization problem: 

WD,opt = argminw^S^RSDWD subject to w^S^s = 1, (1.10) 

■WD 

The optimal low-rank MVDR solution for the above problem is given by 

s" Sd{SpRSd) SjjS sf,Rjj sd 

where Rd = S^RSd denotes the low-rank covariance matrix and sjj = 
S'j^s denotes the low-rank steering vector. The key challenge in the design 
of low-rank STAP algorithms is to find a cost-effective method to compute 
the rank-reduction matrix Sd- 
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1.4.1 Eigenvalue- decomposition-based algorithms 

The eigenvalue-decomposition (EVD)-based beamforming algorithms are 
also known as PC-based algorithms and have been originally reported as 
the eigencanceler method. These PC-based algorithms refer to the beam- 
formers constructed with a subset of the eigenvectors of the interference- 
only covariance matrix associated with the eigenvalues of largest magni- 
tu de. The first appli cation of this method to radar systems was reported 



Haimovich I (|l99ll )]. 



The basic idea of the EVD-based beamformer is to approximate the 
M X M covariance matrix R of the received data as follows: 



D 



(1.12) 



where the M x 1 vector Vd corresponds to the dth eigenvector of R and is 
the dth eigenvalue of R. By assuming that the eigenvalues are obtained in 
decreasing order of magnitude, the EVD-based method approximates R us- 
ing its D dominant eigenvectors. The rank-reduction matrix is constructed 
by using the D dominant eigenvectors as described by 



Sd ^[vi V2 ... Vd] 
The low-rank MVDR solution for the above problem is given by 



Wd = 



s"Sd{S^RSd)-'S^s s^iZti K'-^dV^)^ 



(1.13) 



(1.14) 



The EVD-based low-rank MVDR space-time beamformer described above 
does not take into account the target steering vector s when selecting a 
suitable subspace representation of the interference. Clearly, this low-rank 
space-time beamformer requires the compu tation of an EVD, which h as 
a computational cost that is cubic with M Golub and van Loan ( 20021 )]. 
In order to reduce this computational complexity, a designer can resort 
to sub s pace t racking a l gorith ms which bring the cost down to 0{A'P) 



Yang I (|l995l) ; iBadeau I (|2005[ )] . Another technique associated with EVD- 



based beamforming that can improve the performance of low-rank MVDR 
space-time beamformers is t he method called cross-spectral metric (CSM) 



Goldstein and Reedl (|l997c() ]. The CSM approach chooses the set of D 



eigenvectors for the rank-reduction matrix which optimizes the desired cri- 
terion, namely, the maximization of the SINR, in opposition to the PC 
method which always chooses the dominant eigenvectors. 



February 12, 2013 2:27 



World Scientific Book - 9in x 6in 



ws-book9x6 



Adaptive Space-Time Beamf arming in Radar Systems 



1.4.2 Krylov sub space-has ed algorithms 

The first Krylov methods, namely, the conjugate gradie nt (CG) metho d 



Hestenes and Stiefel I (jl952l )] and the Lanczos algorithm Lanczos I (jl952l )] 



have been originally proposed for solving large systems of linear equations. 
These algorithms used in numerical linear algebra are mathematically iden- 
tical to each other and have been derived for Hermitian and positive def- 
inite system matrices. Other techniques have been repo rted for solving 
these problems and the Arnoldi algorithm Arnoldi I (jl95l[ )] is a computa- 
tionally efficient procedure for arbit rarily invertible syst e m mat rices. The 
multistage Wiener filter (MSW F) IColdstein and Reedl |l997cl)] and the 
auxiliary vector filtering (AVF) Pados and Batalama HJ999)] algorithms 
are based on a multistage decomposition of the linear MMSE estimator. A 
key feature of these methods is that they do not require an EVD and have a 
very good performance. It turns out that Krylov subspace algorithms that 
are used for solving very large and sparse systems of linear equations, are 
highly suitable alternatives for designing low-rank space-time beamforming 
algorithms in radar systems. The basic idea of Krylov subspace algorithms 
is to construct the rank-reduction matrix S^j with the following structure: 



SD=[qRq ■■■ R^'-'q] 



(1.15) 



where q 



Mi\\ 



and 



denotes the Euclidean norm (or the 2-norm) 



of a vector. In order to compute the basis vectors of the Krylov subspace 
(the vectors of So), a designer can either directly employ the expression 
in (jl.lSp or resort to mo re sophisticated approaches such as the Arnoldi 



Arnoldi I (|l95l[ )]. The low-rank MVDR solution for the space- 



iteration 

time beamformer using the Krylov subspace is given by 



Wd 



{S^RSd)-'S^s 



Sd{S^RS 



(1.16) 



An appealing feature of the Krylov subspace algorithms is that the required 
model order D does not scale with the system size. Indeed, when M goes to 
infinity the required D remains a finite and rel atively small value. This re- 
sult was established in [Xiao and Honig I (|2005l )]. Among the disadvantages 



of Krylov subspace methods are the relatively high computational cost of 
constructing So {0{DA'P)), the numerical instability of some implemen- 
tations and the lack of flexibility for imposing constraints on the design of 
the basis vectors. 
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1.4.3 JlO-based algorithms 

The aim of this part is to introduce the reader to low-rank beamforming 
algorithms based on joint iterative optimization (JIO) techniques. The idea 
of these methods is to design the main components of a low-rank space-time 
beamforming scheme via a general optimiz ation approach. The basic ideas 



of JIO t echniques have be en reported i n Ide Lamare and Sampaio-Neto . 
( 2007al): IPa et al\ J2OO8I): Ide Lamare I (|2008l ): Ide Lamare et al. I |20ld ): 



Fa and de Lamare I (|201ll )]. Amongst the advantages of JIO techniques 
are the flexibility to choose the optimisation algorithm and to impose con- 
straints, which provides a significant advantage over eigen-based and Krylov 
subspace methods. One disadvantage that is shared amongst the JIO tech- 
niques, eigen-based and Krylov subspace methods are the complexity asso- 
ciated with the design of the matrix S d ■ For instance, if we are to design a 
beamforming algorithm with a very large M, we still have the problem of 
having to design an M x D rank-reduction matrix S d . 

In the framework of JIO techniques, the design of the matrix S d and the 
beamforming vector wd for a fixed model order D will be dictated by the 
optimization problem and the algorithm chosen to compute the solution. 
To this end, we will focus on a generic S£) — [si S2 ... in which 

the basis vectors s^, d — 1,2, . . . , D will be obtained via an optimization 
algorithm and iterations between the S^i and wd will be performed. The 
JIO method consists of solving the following optimization problem 



['S'_D,opt,i^_D,opt] = arg min w%S"RSdWd, 

x(i) 



subject to w^S^s = 1 



(1.17) 



where it should be remarked that the optimization problem in (jl.l7p is non 
convex, however, the algorithms do not present convergence problems. Nu- 
merical studies with JIO methods indicate that the minima are identical and 
global. Proofs of global convergence have been established wit h different 



versio n s of JIO sche i nes Ide Lamare and Sampaio-Neto I (l2007al): iFa et al. 
( 2008 ): de Lamare ( 20081) : de Lamare et al. I ( 2010l ): Fa and de Lamare 



201 ll )] , which demonstrate that a least-squares (LS) algorithm converges 



to the reduced-rank Wiener filter. 

In order to solve the ab ove problem, we resort to the method of Lagrange 
multipliers HavkinI (|2002[ )] and transform the constrained optimization into 
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an unconstrained one expressed by the Lagrangian 



C{Sd,wd) = wIsIRSdwd + Kwd^dS - 1) 



(1.18) 



where A is a scalar Lagrange multipher. By fixing wo, minimizing p.l^ 
with respect to Sd and solving for A, we obtain 



R 



,"R^^wdsHR-^s' 



(1.19) 



where R — E[r{i)r^ [i)] and Rw — E[w dw^]. By fixing Sdi minimizing 
(|1.18[) with respect to wd and solving for A, we arrive at the expression 



R~ 



s^R- 



(1.20) 



where Rd = E[Sar{i)r" {i)SD] -B[r£,(i)rg(i)], sd = S"s. Note that 
the expressions in (|1.19p and (|1.20p are not closed- form solutions for wjj 
and Sd since (I1.19P is a function of wd and (|1.20l) depends on Sd- Thus, 
it is necessary to iterate (|1.19l) and (|1.20p with initi al values to obtain a 
solution . Unlike the Kr ylov subspace-based me t hods [Goldstein and Reed 
(|l997ah ] and the AVF |Pados and KarvstinosI |200ll )] methods, the JIO 
scheme provides an iterative exchange of information between the low-rank 
beamformer and the rank-reduction matrix and leads to a simpler adaptive 
implementation. The key strategy lies in the joint optimization of the filters. 
The rank D must be set by the designer to ensure appropriate performance 
or can be estimated via another algorithm. In terms of complexity, the JIO 
techniques have a computational cost that is related to the optimization 
algorithm. With recursive LS algorithms the complexity is quadratic with 
M ((0(M^)), whereas the complexity can b e as low as linear with M when 



stocha stic gradient algorithms are adopted jde Lamare and Sampaio-Neto 
|2009h ]. 



1.4.4 JIDF-based algorithms 

This section reviews a low-rank space-time beamforming technique based 
on the joint interpolation, decimation a n d filtering (JID F) concept 



de Lamare and Sampaio-Neto I |2007bll2009l ) : iFa et al. I |2010l )]. The JIDF 



approach allows a designer to compute the parameters of the rank- reduction 
matrix and the low-rank space-time beamformer with a low complexity. 
The motivation for designing a rank-reduction matrix based on interpola- 
tion and decimation comes from two observations. The first is that rank 
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reduction can be performed by constructing new samples with interpola- 
tors and eliminating (decimating) samples that are not useful in the STAP 
design. The second comes from the structure of the rank-reduction matrix, 
whose columns are a set of vectors formed by the interpolators and the 
decimators. 

In the JIDF scheme, the number of elements for adaptive processing 
is substantially reduced, resulting in considerable computational savings 
and very fast convergence performance for the radar applications. The 
M X 1 received vector r{i) is processed by a multiple processing branch 
(MPB) scheme with B branches, where each spatio-temporal processing 
branch contains an interpolator, a decimation unit and a low-rank space- 
time beamformer. In the 6-th branch, the received vector r(i) is filtered by 
the interpolator Vb = [vb,o Vb,i ■ . ■ Vb^i-i]'^ with / coefficients, resulting in 
an interpolated received vector r(,(i) with M samples, which is expressed 
by 

nil) = r(z), (1.21) 
where the M x M Toeplitz convolution matrix is given by 

... 



Vb,0 
VbJ-l 








Vb,I-l 






"6,0 



(1.22) 



The vector rb{i) can be expressed in an alternative way that is useful for 
the design of the JIDF scheme and is described by 

nil) = V^r{i) = Ua{i)vb, (1.23) 

where the M x I matrix Sto{i) with the samples of r{i) has a Hankel struc- 
ture and is described by 

ro(?:) ri{i) ...ri_i(i)' 
ri{i) r2{i) ... n(«) 



5Ro(«) 



rM-2{i) rM-i{i) 
.rM-i{i) 



(1.24) 
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(1.25) 



The dimensionality reduction is performed by a decimation unit with DxM 
decimation matrices Db that transforms r/(i) into D x 1 vectors rD,b{i) 
with b = 1, . . . ,B, where D — AI/L is the rank of the resuhing system of 
equations that will be generated and L is the decimation factor. The 13 x 1 
vector r]j^h{i) for branch b is expressed by 

rD,b - ^g,fcr(z) = DD,bV^r{i) 

where SD,b is the rank-reduction matrix and the vector r£),b(z) for branch 
b is used in the minimization of the output power for branch b. The output 
at the end of the JIDF scheme y{i) is selected according to 

yii) = VbAi) when 6^ = argmin ly^p (1-26) 

b 

where i? is a parameter to be set by the designer. For the computation of 
the parameters of the JIDF scheme, it is fundamental to express the output 
yb{i) as a function of the interpolator Vh, the decimation matrix Dd^i, and 
the low-rank space-time beamformer wjj^i, as follows: 

yb{i) = w^j,S"^br{i) 

= w",^DD,b^oii)v, 

where the expression (|1.27p indicates that the dimensionality reduction car- 
ried out by the JIDF scheme depends on finding appropriate ■»(,, DD,b and 
WD,b- Unlike the previously discussed low-rank beamforming techniques, 
the JIDF is able to substantially reduce the cost of the rank-reduction ma- 
trix. 

The parameters of the JIDF scheme that perform low-rank space-time 
MVDR beamforming can be computed by solving the following optimiza- 
tion problem 

[iOD,opt, ■"opt, -C>D,bJ = arg min w" ^E[DD,b^o{i)vv"^" {i)D'j^ ^]wd 
subject to w^j^Do.bSoVb — 1, 

(1.28) 

where So is M x I steering matrix with a Hankel structure, which has the 
same form as 5io(i) and is given by 

so(«) si{i) . . . si_i(i) 
Sl(j) S2(«) ••■ si(j) 



(1.27) 



So{^) 



Sj\/-2(j) SAf-l(j) ... 
SM-l(j) ... 



(1.29) 
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The constrained optimization in (|1.28|) can be transformed into an uncon- 
strained optimization problem by using the method of Lagrange multipliers, 
which results in 

£{wD,b,Vb,DD.b) = w'^E[DDMi)vv"U^ {i)D'^ .jwD+Kw'^DD^bSoV-l), 

(1.30) 

where A is a Lagrange multiplier. 

The strategy to compute the parameters of the low-rank space-time 
beamformer based on the JIDF scheme is to minimize the cost function 
with respect to a set of parameters and fix the remaining parameters. By 
minimizing (|1.30p with respect to v^, we obtain 

R-l„ 

v,b''v,b 
''v,b-^v,b''v,b 

where Rv,b — £'[r^,br^j] is the I x I autocorrelation matrix, r^^j, = 
-D^ ^R^WBfi, and s^.b — i^S^wdj, is the Ix low-rank steering vector. 
By minimizing (jl.30p with respect to wo.b^ we have 

Rw,b^it',b 

WD,b= 1 , (1-32) 

^w,b^w,b^w,b 

where Rw,b — Elr^^fwi^ '^^ D x D autocorrelation matrix, r^^b = 
DD,bRoVb, and Sw^b — Do.bSoVb is the D x I low-rank steering vector. In 
order to compute Vb and Wjj bi a designer needs to iterate them for each 
processing branch b. 

The decimation matrix b is selected to minimize the square of the 
output of the beamformer yb{i) obtained for all the B branches 

Do.b^Do.bM when 6s arg min \yb{i)\'^ , (1.33) 

l<b<B 

The design of the decimation matrix D u b imposes constraints on the values 
of the elements of the matrix such that they only take the value zero or 
one. Since the optimal approach for the design of Djj b corresponds to an 
exhaustive search, we consider a suboptimal technique that employs pre- 
stored patterns. The decimation scheme employs a structure formed in the 
following way 

So.b = [</'64 0b,2 Am], (1-34) 
where <pb d an M x 1 vector composed of a single one and zeros elsewhere 
as described by 

ct>b^^ = 1, 0^_^], (1.35) 
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where Zt^d is the number of zeros before the only element equal to one. We 
set the value of zt^d in a deterministic way which can be expressed as 

zt,d = ^^(d-l) + ib-l). (1.36) 

It is necessary to iterate (|1.3ip . (I1.32[) and (ll.33|) in an alternated form 
(one followed by the other) with an initial value to obtain a solution. The 
expectations can be estimated either via time averages or by instantaneous 
estimates and with the help of adaptive algorithms. 



1.5 Sparsity- Aware Beamforming Algorithms 

This section considers space-time beamforming algorithms that exploit the 
sparsity encountered in the data processed by radar systems. In particu- 
lar, the motivation for exploiting the sparsity of data vectors observed by 
radar systems is given and a brief discussion on the suitability of sparsity- 
aware algorithms for radar applications is provided. A general approach 
to design space-time beamforming algorithms based on the Zi-norm regu- 
larization is described. The main principle is to employ a reduced number 
of weights to suppress the clutter and the jamming signals encountered in 
radar applications. 

Recently, motivated by compressive sensing (CS) techniques used in 
radar, several authors have consi dered CS ideas for mov i ng target indicatio n 
(MTI) and STAP problems [Mar ia and Fuchs I |2006l )]- |Selesnick I |2O10l )]. 



The core notion in CS is to regularize a linear inverse p roblem by including 



prior knowledge that the signal of interest is sparse Parker and Potter 



(|2010f )]. These works on space-time beamforming techniques based on CS 
rely on the recovery of the clutter power in angle-Doppler plane, which is 
usually carried out via convex optimization tools. However, these methods 
are based on linear programming and have a quite high computational 
complexity {0{K^)), where K is the dimension of the angle-Dopplcr plane. 
In this section, we describe the concept of a sparsity-aware STAP (SA- 
STAP) algorithm that can improve the detection capability using a small 
number of snapshot. To overcome the high complexity of the CS-STAP 
type algorithm, we design the STAP algorithm with another strategy, by 
imposing the sparse regularization to the minimum variance (MV) cost 
function. Since the interference variance has often a low-rank property, we 
assume that a number of samples of the data cube are not meaningful for 
processing and the optimal STAP beamformer is sparse, or nearly sparse. 
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Then, we exploit this feature by using a Zi-norm regularization. With this 
motivation, the STAP algorithm design becomes a mixed /i-norm and I2- 
norm optimization problem. 

The conventional space-time beamforming algorithms do not exploit 
the sparsity of the received signals. In this exposition, it is assumed that 
a number of samples of the data cube are not meaningful for processing 
and a reduced number of active weights of the space-time beamformer can 
effectively suppress the clutter and the jamming signals. Specifically, a 
sparse regularization is imposed to the space-time MVDR beamforming 
design. Thus, the space-time beamformer design can be described as the 
following optimization problem 

Wont — avgimnw^Rw 

- (1.37) 
subject to w s = 1 and \ \w\\i = 0, 

where the objective of the Zi-norm regula rization is to force th e comp onents 



of the space-time beamformer w to zero jAngelosante et al. I pOlOf l]. This 
problem can be solved using the method of Lagrange multipliers, which 
results in the following unconstrained cost function 

C{w,a,X) = w" Rw + a{w" s - 1) + X{\\w\\i), (1.38) 

The unconstrained cost function above is convex, however, it is non- 
differentiable which makes it difficult for one to use the method of La- 
grange Multipliers directly and obtain an expression for the space-time 
beamformer. To this end, the following approximation to the regulariza- 
tion term is employed 

\\w\\i^w"Aw, (1.39) 

where 

A - diag f . . . I \—) , (1.40) 

where e is a small positive constant. Simultaneously, we assume that the 
partial derivative of Aw with respect to w* is given by 

-g^^Au,. (1.41) 

With the development above, an approximation to the unconstrained cost 
function can be employed as described by 

C{w, a, A) « w"Rw + a{w"s - 1) + Xw" Aw, (1.42) 
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By computing the gradient terms with respect to w* and a and equat- 
ing them to zero, we obtain the following expression for the space-time 
beamformer 

w = ^— /. , , . 1.43 

Comparing (jl.43l) with the conventional optimal space-time beamformer 
in (|1.8|) . we find that there is an additional term A A in the inverse of the 
interference covariance matrix ii, which is due to the Zi-norm regulariza- 
tion. The term A is a positive scalar which provides a trade-off between 
the sparsity and the output interference power. The larger the chosen A, 
the more components are shrunk to zero [Zibulevskv and Elad ] lioio)]. It 
should also be remarked that the expression for the beamformer in ()1.43p 
is not a closed-form solution since A is a function of w. Thus it is neces- 
sary to develop an iterative procedure to compute the parameters of the 
space-time beamformer. 



1.6 Knowledge- Aided Beamforming Algorithms 



Although STAP techniques are considered efficient tools for detection of 
slow targets b y airborne radar systems in strong clutter environments 
Klemml(|2003 l]. due to the very large number of degrees of freedom (DoFs) 



conventional space-time beamformers have a slow convergence and re- 
quire about twice the DoFs of the independent and identically distributed 
(IIP) training snap shots to yield an average performance loss of roughly 
3dB Ward I ([l99J)]. In real scenarios, it is hard to obtain so many IID 
training snaps h ots, especially in heterogeneou s envir onments. Low-rank 
Guerci ( "2000[l1- de Lamare and Sampaio-Neto ( 20091 )] and sparsity-aware 



Maria and Fuchs 



( 20061 1]- [Yang et al. I (|201lh ] methods have been consid- 



ered to counteract the slow convergence of the conventional space-time 
beamformers. Nevertheless, there are other alternatives to improve the 
training of STAP algorithms and improve their performance. These other 
methods can also be combined with the techniques previously discussed. 
Recently developed knowledge-aided (KA) STAP algorithms have received 
a growing interest and become a key conce p t for the next ge neration of 



a growing interest ana pecome a Jtey conce p t tor tne next ge neration ot 
adaptive radar systems Wicks et al. I ( 20061 )]- Fa et al. ( 2O10l )]. The core 
idea of KA-STAP is to incorporate prior knowledge, provided by digital 
elevation maps, land cover databases, road maps, the Global Position- 
ing System (GPS), previous scanning data and other known features, to 
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compute estimates of the clutter covariance matrix with high accuracy 



Melvin and Guerci I (j2006t )]. Prior work on KA-STAP algorithms include 



the exploitati on of prior knowledge of the clutter ridge to form the STAP 
filter weights Melvin and Showman ( 2006I )]. use of prior knowledge about 



the terrain [Capraro et al. I (|2006[ )] and prior knowle dge about th e covar l 



ance matr i x of t he clutter and the jamming signals [Blunt et al. ( 20061) ]- 



Fa et al. 



In this section, we discuss a strategy to mitigate the deleterious effects 
of the heterogeneity in the secondary data, which makes use of a priori 
knowledge of the clutter covari ance matrix and has r ecently gained si gnifi- 



cant attention in the literature (Wicks et al7\ (|2006[ )]- |Fa et al. I |2010l )]. In 



KA-STAP techniques, there are two basic tasks that need to be addressed: 
the first one is how to obtain prior knowledge from the terrain knowledge 
of the clutter and how to estimate the real interference covariance matrix 



with the prior knowledge [Wicks et al. I (|2006l )]- [Capraro et al. I (|2006l )] and 
the second is how to apply the covariance matrix estimates in the design 



of the space-time beamforming algorithm Blunt et al. I (|2006r )]- Fa et al. 
(|2010l )]. We first review how a designer can obtain prior knowledge of the 



clutter and employ this knowledge to build a known covariance matrix i?o ■ 
Then, we present a method to combine this prior knowledge with com- 
monly used estimation techniques to compute the covariance matrix of the 
received vector r{i), resulting in a combined covariance matrix estimate 
Rc for use in the space-time beamformer that is more accurate and has an 
enhanced performance. 

The optimal space-time beamformer employs the following expression 
to compute its parameters 



R 



s"R 's' 



(1.44) 



where an estimate of the covariance matrix is typically obtained by 



(1.45) 



k=l 



where r{k) is taken from secondary data. The e stimate R can be sufficientl y 
accurate when K is at least twice as great as M Brennan and Reed I (|l973r )] 
and the training samples are assumed i.i.d. However, it is by now well 
understood that the clutter environments are often heterogeneous and this 
leads to performance degradation on space-time beamforming. KA-STAP 
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techni ques can significantly help to combat the heterogeneity [Stoica et al. 



torn ]. 

With KA techniques the covariance matrix is estimated by combin- 
ing an initial guess of the covariance matrix Ro derived from the digital 
terrain database or the data probed by radar in previous scans, and the 
sample average covariance matrix estimate in the present scan R so that 

Rc = aRo + (1 - a)R, (1.46) 

where < a < 1 . Alternatively, this principle can be applied to the inverse 
of the covariance matrix estimate 

R-' ^ fjR-' + {1 - rj)R~\ (1.47) 

where < < 1. 

In order to compute the parameter 77, we need to consider the optimiza- 
tion problem 

Vopt = argmin Rw, (1-48) 

where we use the relation 

w{i) = 1TW0 + {1 - v)w{i), (1.49) 
where w{i) — R^s, Wo — Rq^^ and w — R s. We can obtain the 



optimal value for rj by ec 
to zero, which results in 



uating the gr a dient of the cost function in (jl.48p 



Stoica et al. I (|2008[ )] 



^{s"{R ^ R-^)RR ^s) 



Vopt = ^ / _i ° \ (1-50) 

s"{R-^-R )R{R-^ R )s 

Since R above is unknown, we have to estimate it in real time using either 
time averages or adaptive algorithms. 



1.7 Simulations 

In this section, we assess the performance of the space-time beamforming 
algorithms discussed in this chapter using simulated radar data. Specifi- 
cally, we consider the optimal MVDR space-time beamforming algorithm 
that assumes perfect knowledge of the covariance matrix of the received 
data, and the MVDR space-time beamformer using the sample matrix in- 
version (SMI-MVDR). The low-rank space-time beamforming algorithms 
using EVD (LR-EVD), the Krylov subspace approach (LR-Krylov), the 



February 12, 2013 2:27 



World Scientific Book - 9in x 6in 



ws-book9x6 



20 



My Book Title 



JIO (LR-JIO) and the JIDF (LR-JIDF) algorithms are also considered with 
a rank equal to D. We also consider the sparsity-aware (SA-MVDR) and 
the knowledged-aided (KA-MVDR) space-time beamforming algorithms. 
All the analyzed algorithms estimate the statistical quantities via time- 
averages in a similar way to a least-squares method. The parameters of the 
simulated radar platform are shown in Table 11.11 For all simulations, we 
assume the presence of a mixture of two broadband jammers at —45° and 
60° with jammer-to- noise ratio (JNR) equal to 40 dB. The clutter-to-noise- 
ratio (CNR) is fixed at 40 dB. All the results presented are averages over 
1000 independent Monte-Carlo runs. 



Parameter Value 

Antenna array Sideway-looking array (SLA) 
Carrier frequency (/c) 450 MHz 

Transmit pattern Uniform 



In the first experiment, we assess the SINR performance of the different 
space-time beamforming algorithms as shown in Fig. 11.31 The algorithms 
are simulated over 800 snapshots and the SNR is set to 10 dB. The re- 
sults show that the LR-JIDF algorithm achieves the best results, followed 
by the LR-JIO, the KA-MVDR, the SA-MVDR, the LR-Krylov, the LR- 
EIG and the SMI MVDR algorithms. The curves indicate that the use of 
low-rank algorithms is highly beneficial to the performance of space-time 
beamforming algorithms in radar systems. In particular, the LR-JIDF and 
LR-JIO algorithms have a very fast convergence performance. It should 
also be remarked that the SA-MVDR and KA-MVDR algorithms obtain a 
performance that is significantly better than the conventional SMI MVDR 
algorithm. Since the SA-MVDR and KA-MVDR techniques are modifica- 
tions of the SMI MVDR techniques exploiting sparsity and prior knowledge 
about the covariance matrix, respectively, it is interesting to note that by 
exploiting these properties it is possible to significantly improve the perfor- 
mance of beamforming algorithms. 




75 m/s 
9000 m 
40 dB 



300 Hz 



8 
8 
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Snapshots 

Fig. 1.3 SINR Performance against the snapshots. Parameters; B = 8, I = 8 and 
D = G, M = 64, R[0] = 0.017, e = 0.1. 



We evaluate the SINR performance against the target Doppler frequency 
at the main bean look angle for our proposed algorithms and other exist- 
ing algorithms, which are illustrated in Fig. 11.41 The potential Doppler 
frequency space form —100 to 100 Hz is examined and 100 snapshots are 
used to train the beamformers. The plots show that the analyzed algo- 
rithms converge and approach the optimum in a short time, and form a 
deep null to cancel the main beam clutter. Again, the LR-JIDF algorithm 
outperforms the other analyzed algorithms. 

In the third example, the probability of detection Pd versus SNR per- 
formance is presented for all schemes using 200 snapshots as the training 
data as shown in Fig. 11.51 The false alarm rate PFA is set to 10~^ and we 
suppose the target injected in the boresight (0°) with Doppler frequency 
lOOHz. The figure illustrates that the analyzed algorithms provide sub- 
optimal detection performance using short support data. Note that for 
Pd — 0.9 (90% percent), the LR-JIDF and LR-JIO schemes are within less 
than 1 dB from the performance of the optimal MVDR algorithm. The 
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30 




Target Doppler Frequency ^ (Hz) 



Fig. 1.4 SINR Performance against the target Doppler frequency. Parameters: B = S, 
/ = 8 and D = 6, M = 64, R[0] = 0.01/, e = 0.1. 



remaining techniques exhibit increasing performance losses as compared to 
the optimal MVDR algorithm and it should be noted that the conventional 
SMI MVDR method has a performance degradation of up to 5 dB for the 
same performance measured in terms of Pp. This suggest that the ap- 
plication of more sophistical space-time beamforming algorithms is key to 
achieving an improved performance. 

At last, we detail the computational complexity in terms of multipli- 
cations of the analyzed schemes in Fig. 11.61 The curves show that the 
computational complexity of the LR-JIDF and LR-Krylov algorithms is 
significantly lower than the remaining algorithms. Indeed, there is a sig- 
nificant computational advantage obtained by using the LR-JIDF and LR- 
Krylov algorithms and this advantage becomes more pronounced as M is 
increased. The other analyzed algorithms have a complexity that scales as 
a cubic function of M. This high complexity can be mitigated by the use of 
adaptive algorithms, which can reduce the computational cost by at least 
an order of magnitude. 
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Fig. 1.5 Probability of detection against SNR. Parameters: B 
M = 64, R\0] = O.Oli", e = 0.1. 



and D = 6, 



1.8 Concluding Remarks 



This chapter has presented recent advances on space-time beamforming 
algorithms for phased-array radar systems and investigated their perfor- 
mance via computer simulations. Specifically, we have considered in detail 
some of the most successful space-time beamforming algorithms that ex- 
ploit low-rank and sparsity properties as well as the use of prior-knowledge 
to improve the performance of STAP algorithms. The results of our studies 
suggest that the low-rank algorithms have a substantial performance advan- 
tage over conventional MVDR space-time beamforming algorithms. Fur- 
thermore, the use of sparsity- aware and knowledge-aided strategies is also 
able to improve the performance of space-beamforming algorithms, and can 
be combined with low-rank schemes. These beamforming algorithms can 
be also applied to MIMO radar systems, mono-static and bi-static radar 
systems and other sensing applications such as sonar systems. 
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Fig. 1.6 Computational complexity in terms of multiplications of the analyzed space- 
time beamforming algorithms. Parameters: S = 8, / = 8 and _D = 6. 
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